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Abstract 

The 4 He total photoabsorption cross section is calculated with the realistic nucleon-nucleon 
potential Argonne V18 and the three-nucleon force (3NF) Urbana IX. Final state interaction is 
included rigorously via the Lorentz Integral Transform method. A rather pronounced giant res¬ 
onance with peak cross sections of 3 (3.2) mb is obtained with (without) 3NF. Above 50 MeV 
strong 3NF effects, up to 35%, are present. Good agreement with experiment is found close to 
threshold. A comparison in the giant resonance region is inconclusive, since present data do not 
show a unique picture. 

PACS numbers: 21.45.+v, 21.30.Fe, 25.20.Dc, 24.30.Cz 
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In the last three decades there has been a continous interest in 4 He photodisintegration, 


both in theory and in experiment (see 


IS 


and references therein). The a-particle is 


drawing such a great attention because it has some typical features of heavier systems (e.g. 
binding energy per nucleon), which make it an important link between the classical few- 
body systems, i.e. deuteron, triton and 3 He, and more complex nuclei. For example in 
4 He one can study the possible emergence of collective phenomena typical of complex nuclei 
like the giant resonance. Furthermore 4 He is the ideal testing ground for microscopic two- 
and three-body forces, which are fitted in the two- and three-body systems. At present the 
3NF is not yet well determined, thus it is essential to search for observables where it plays 
an important role. Because of gauge invariance, in electromagnetic processes nuclear forces 
manifest themselves also as exchange currents, which have turned out to be very important 


le three- 
it was 


1C IJ 

I it 


in photonuclear reactions and hence 3NF effects might become significant. For t 
nucleon systems photonuclear processes have already been studied [ 111 , 0 ], e.g. in 
found that the 3NF leads to an almost 10% reduction of the electric dipole peak and up to 
15% enhancement at higher energy. One expects that the 3NF is of considerably greater 
relevance in the four-body system, since there one has six nucleon pairs and four triplets, 
compared to three pairs and just one triplet in the three-nucleon systems. 

Unfortunately, the current theoretical and experimental situation of the 4 He photodisin¬ 
tegration is not sufficiently settled. Most of the experimental work has concentrated on the 
two-body break-up channels 1 He( 7 , n) 3 He and 4 He( 7 ,p) 3 H in the giant resonance region, but 
still today there is large disagreement in the peak. In fact in two very recent ( 7 , n) experi¬ 
ments 0, l^l one finds differences of a factor of two. The 4 He( 7 ) reaction represents a very 
challenging theoretical problem as well, since the full four-body continuum dynamics and 
all possible fragmentations have to be considered. Calculations with realistic nuclear forces 
have not yet been carried out. Even investigations with a realistic nucleon-nucleon (NN) 
interaction and without a 3NF are still missing. Calculations exist only for semi-realistic 


NN potentials. In Refs. 


[I II B 


it has been shown that such models lead to pronounced 


0 


peak cross sections, in rather good agreement with the experimental data of 1 ] and much 


different from what was calculated earlier 

It is evident that the experimental and theoretical situations are very unsatisfactory. 
With the present work we make an important step forward on the theory side performing 
a calculation of the a-particle total photoabsorption cross section with a realistic nuclear 
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force. To this end we solve the four-body problem taking as nuclear interaction the realistic 
Argonne V18 (AV18) NN potential 111 and the Urbana IX (UIX) 3NF [u|. We reduce the 
continuum state problem to a bound-state like problem via the Lorentz Integral Transform 
(LIT) method jlOj], thus taking into account the full final state interaction rigorously. The 
LIT bound-state like problem, as well as that of the 4 He ground state, is solved via expan¬ 
sions in hyperspherical harmonics (HH) using the powerful effective interaction HH (EIHH) 


ia 


14]. The reliability of the EIHH technique in electromagnetic processes with 


approach 

the AV18 and UIX potentials has been shown in a recent 3 H(y) calculation |l5|. 

We calculate the total photoabsorption cross section <r 7 (u;) in the unretarded dipole ap¬ 
proximation. In this way the dominant part of the exchange current contribution on the 
total cross section is taken into account via the Siegert theorem. In the classical few-body 
systems it has been found Jj[ [u| that this is an excellent approximation, particularly for 
photon energies u below 50 MeV. For example for the triton <r 7 (u;) the contributions of 
retardation and all other multipoles lead to a cross section enhancement of less than 1% for 
uj < 40 MeV and of 5 (5)%, 16 (18)%, and 26 (33)% with AV18 (AV18+UIX) at uj — 60, 
100, and 140 MeV, respectively [4]. 

The total photoabsorption cross section is given by 


cr 7 (u;) = 4:TT 2 aujR(uj ), 


( 1 ) 


where a is the fine structure constant and 


R(u>) = / d^f ($ f \D z \V 0 ) 5(E f -E 0 -u) 


( 2 ) 


is the response function in the unretarded dipole approximation with D z = -^p. The 

wave functions of the ground and final states are denoted by | v L 0 /j) and the energies by 
E 0 /f, respectively. The operators rf and z[ are the third components of the i-th nucleon 


isospin and center of mass frame position. In the LIT met 
inversion of an integral transform with a Lorentzian kernel 

R(u) 


lod one obtains R(uj) from the 
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= (TO 


L(<jr, oy) — / du- 2 

J (uj- a R ) 2 + erf 

where T is the unique solution of an inhomogeneous “Schrodinger-like” equation 


( 3 ) 


(H - E 0 - a R + hT/)|T) = D z |T 0 ) 


( 4 ) 
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with bound-state like asymptotic boundary conditions. 

The EIHH expansions of To and T are performed with the full HH set up to maximal 
values of the HH grand-angular momentum quantum number K (K < A// for To, K < K m 
for T). The convergence of binding energy and matter radius are presented in Table I. 
Our final binding energy results of 24.27 MeV (AV18) and 28.42 MeV (AV18+URBIX) 
agree quite well with other calculations (AV18: 24.25 171, 24.22 24.21 JhJ MeV); 

AV18+URBIX: 28.34 [20], 28.50 171, 28.46 Q MeV). 

The EIHH convergence of the transform £ is excellent for the AV18 potential. Here we 
discuss in detail only the case AV18+UIX, where the convergence is quite good, but not 
at such an excellent level. The reason is that in our present EIHH calculation an effective 
interaction is constructed only for the NN potential, while the 3NF is taken into account as 
bare interaction. In Fig. 1 we show results for the transform £ obtained with various K m 
and K Since we take K Q rn = K rn — 1 the corresponding transform can be denoted by L Km . 
One sees that there is a very good convergence beyond the peak. The figure also shows that 
the peak height is very well established, but that the peak position is not yet completely 
converged. In fact with increasing K m the peak is slightly shifted towards lower a R . This is 
illustrated better in Fig. 2, where we show the relative differences A Km = L Kml9 /Li 9 with 
L a ^ = L a — Lp (the chosen OR-range starts at the 4 He( 7 ) break-up threshold). One again 
notes the very good convergence for cr R > 30 MeV with almost identical results from £13 
to L w . Altogether we consider our result for a R > 30 MeV as completely sufficient. On 
the other hand it is obvious that convergence is not entirely reached for lower ctr. Here 
we should mention that for L\ 9 there are already 364000 states in the HH expansion. A 
further increase is beyond our present computational capabilities. On the other hand, as a 
closer inspection of Fig. 2 shows, the convergence proceeds with a rather regular pattern: 
(i) £ 13,11 — T 11,9 and £ 17,15 ~ £ 15,13 and (ii) £ 19,17 — £ 17 , 15 / 1.5 ~ £i3,n/(1.5) 2 . Due to 
this regular pattern it is possible to obtain an extrapolated asymptotic result. We use the 
following Pade approximation 


Lk =L 


A' m _ s + L Km _ 4yKj 


,- 8/(1 


j K K —4 

Jv m 5 Jv m ^ 


Lk 


4 ,Km-8 


)• 
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In Fig. 3 we illustrate the results for obtained from the inversions of the correspond¬ 
ing transforms £/<- m . Due to the Lorentz kernel the cross section presents the same features 
as £ itself: established peak height with a value very close to 3 mb, and not yet completely 
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convergent peak position. In Fig. 3 we also show cr“ 17 and cr“ 19 resulting from the inversion 
of L ^ and respectively. Unquestionably, the extrapolated have a lower numerical 

quality than the calculated Lx m and consequently we do not find the same stability for the 
corresponding inversions (for details on inversion see 21]). Therefore we use an additional 
constraint in the inversion by fixing the peak cross section to the already converged value of 
3 mb. In Fig. 3 it is evident that cx“ 17 and cr“ 19 are very similar, hence establishing a very 
good approximation for the asymptotic <r 7 . One also notices that compared to cr 7i ig they 
show a shift of the peak position by about 1 MeV towards lower energy. 

In Fig. 4a we show our final results for cr 7 . Due to the 3NF one observes a reduction of 
the peak height by about 6% and a shift of the peak position by about 1 MeV towards higher 
energy. Very large effects of the 3NF are found above 50 MeV with an enhancement of er 7 
by e.g. 18, 25, and 35% at uj — 60, 100, and 140 MeV, respectively (3NF effect could change 
somewhat if all multipole contributions are considered, see also discussion before Eq. (1)). 
The comparison between the present results and those obtained with semi-realistic potential 
models p, 7] is similar as in the three-nucleon photo disintegration : semi-realistic models 
lead to rather realistic results in the giant resonance region overestimating the peak by about 
10-15% and giving quite a correct result for the the peak position; however, at higher u the 
cross section is strongly underestimated (in the present case by a factor of three at pion 
threshold). 

It is very interesting to compare the 3NF effects on cr 7 to those found in the three-body 
nuclei |2i, 4J. Surprisingly, the reduction of the 4 He peak height is smaller than in the three- 
nucleon case. For 3 H/ 3 He its size is similar to the binding energy increase (10%), whereas 
for 4 He the 3NF increases the binding energy by 17%, but reduces the peak by only 6% and 
thus cannot be interpreted as a simple binding effect. Also at higher u there are important 
differences. The enhancement of the 4 He cross section due to the 3NF is significantly higher, 
namely about two times larger than for the three-body case. Interestingly this reflects the 
above mentioned different ratios between triplets and pairs in three- and four-body systems. 

In Fig. 4b we compare our results to experimental data. For the AV18+UIX case up¬ 
per/lower bounds are included to account for possible errors in the extrapolation (|SJ). The 
bounds are obtained enhancing/reducing the difference cr“ 19 — ct 7i i 9 by 50%; we believe that 
this is a rather safe error estimate. The data in Fig. 4b have to be interpreted with some 


care: l m 


22] the peak cross section is determined from Compton scattering via dispersion 


5 





relations, (ii) the clashed area corresponds to the sum of cross sections for ( 7 , n) from 


and ( 7 ,p) 3 H from | 24j as already shown in jf|, (iii) the data from the above mentioned recent 
( 7 , n ) experiment 1 ] are included only up to about the three-body break-up threshold, where 


one can rather safely assume that <r 7 ~ 2 cr( 7 , n) (see also | 8 (), (iv) in | 2 | all open channels are 
considered. One sees that the various data are quite different exhibiting maximal deviations 


of a 
0 


rout a factor of two. The theoretical cr 7 agrees quite well with the low-energy data of 


241 ]. In the peak region, however, the situation is very unclear. There is a rather go_od 


agreement between the theoretical <j 7 and the data of 


1] and 221, while those of 23, 24;] 


are noticeably lower. Very large discrepancies are found in comparison to the recent data of 


l a. 


Shirna et al. | 2 |. It is evident that the experimental situation is rather unsatisfactory and 
further improvement is urgently needed. 

We summarize our work as follows. We have performed an ab initio calculation of the 4 He 
total photoabsorption cross section using a realistic nuclear force (AV18 NN potential plus 
the UIX-3NF). The full interaction is consistently taken into account both in the ground state 
and in the continuum via the LIT method. For the solutions of the differential equations 
we use expansions in hyperspherical harmonics via the EIHH approach. Our results show a 
rather pronounced giant dipole peak. The 3NF reduces the peak height by only about 6 %, 
less than expected considering its large effect of almost 20% on the 4 He binding energy and 
its different role in the three-nucleon system. Beyond the giant dipole resonance 3NF effects 
become much larger. With growing u their importance increases and at pion threshold one 
has an enhancement of 35%, about twice the effect one finds in 3 H/ 3 He photoabsorption. 
Close to threshold the theoretical cross section agrees quite well with experimental data. 
In the giant resonance region, where there is no established experimental cross section, 


0 - 


our results are in good agreement with the data of [ 1 ] and 


22 j, while we find a strong 


disagreement with the data of 

In conclusion we would like to emphasize that it is very important to understand whether 
a nuclear force model, which is constructed in the two- and three-nucleon systems, is suf¬ 
ficient to explain the four-nucleon photodisintegration. To this end further experimental 
investigations are mandatory. 
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FIG. 4: Total 4 He photoabsorption cross section: (a) cr 7 (AV18) and <t“ 19 (AV18+UIX), (b) as 
(a) but also included upper/lower bounds and various experimental data (see text), area between 


dotted lines P, 24] , dotted box 


22], squares 


Q|, and circles 
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TABLE I: Convergence of HH expansion for 4 He binding energy E b [MeV] and root mean square 
matter radius (r 2 }5 [fm] with AV18 and AV18+UIX potentials. 


K° 

m 

AV18 

E b 

(r 2 ) 2 

AV18+UIX 

E b 

(r 2 )5 

6 

25.312 

1.506 

26.23 

1.456 

8 

25.000 

1.509 

27.63 

1.428 

10 

24.443 

1.520 

27.861 

1.428 

12 

24.492 

1.518 

28.261 

1.427 

14 

24.350 

1.518 

28.324 

1.428 

16 

24.315 

1.518 

28.397 

1.430 

18 

24.273 

1.518 

28.396 

1.431 

20 

24.268 

1.518 

28.418 

1.432 
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